Design and analysis for the SPICE parameters of waveform-selective metasurfaces varying with the incident pulse width at a constant oscillation frequency

In this study, we numerically demonstrate how the response of recently reported circuit-based metasurfaces is characterized by their circuit parameters. These metasurfaces, which include a set of four diodes as a full wave rectifier, are capable of sensing different waves even at the same frequency in response to the incident waveform, or more specifically the pulse width. This study reveals the relationship between the electromagnetic response of such waveform-selective metasurfaces and the SPICE parameters of the diodes used. In particular, we draw conclusions about how the SPICE parameters are related to (1) the high-frequency operation, (2) input power requirement and (3) dynamic range of waveform-selective metasurfaces with supporting simulation results. First, we show that reducing a parasitic capacitive component of the diodes is important for realization of the waveform-selective metasurfaces in a higher frequency regime. Second, we report that the operating power level is closely related to the saturation current and the breakdown voltage of the diodes. Moreover, the operating power range is found to be broadened by introducing an additional resistor into the inside of the diode bridge. Our study is expected to provide design guidelines for circuit-based waveform-selective metasurfaces to select/fabricate optimal diodes and enhance the waveform-selective performance at the target frequency and power level. Our results are usefully exploited to ensure the selectivity based on the pulse duration of the incident wave in a range of potential applications including electromagnetic interference, wireless power transfer, antenna design, wireless communications, and sensing.

Artificially engineered structures referred to as metamaterials and metasurfaces enable ready tailoring of electromagnetic fields at will [1][2][3][4][5] . Compared to conventional materials that respond to an incident wave in accordance with the reaction of composite molecules, the response of metasurfaces is characterized by composite unit cells that are much larger than molecules but still at the subwavelength scale and can thus be arbitrarily customized to achieve a range of electromagnetic properties including negative refraction 6 and a large impedance surface 3 . In addition, these metasurfaces with exotic properties can be exploited to develop applied devices or systems such as invisibility cloaks 7,8 , diffraction-limit-breaking lenses 9,10 , perfect absorbers [11][12][13][14] , antennas 15,16 , analogue computation systems [17][18][19] and beamforming systems 4,5,[20][21][22] as intelligent reflecting surfaces (IRSs) [23][24][25] .
In particular, these artificially constructed materials attain a higher performance level by including nonlinearity, which offers capabilities that are not available with a simple combination of linear media or structures [26][27][28][29] . For instance, nonlinear metasurface-based absorbers dissipate the energy of a strong electromagnetic wave to Scientific Reports | (2023) 13:7202 | https://doi.org/10.1038/s41598-023-34112-z www.nature.com/scientificreports/ prevent electromagnetic interference issues while permitting propagation of a small signal for antenna communication even at the same frequency [30][31][32][33] . Also, nonlinear metasurfaces can be designed as intensity-dependent surfaces allowing various functionalities including spatial wave control [34][35][36] , switchable transmission 37 and digitally defined beam deflection 38 . Additionally, the electromagnetic response of nonlinear metasurfaces or IRSs is altered in accordance with surrounding electromagnetic fields or incoming communication signals, which is achieved by including nonlinear circuit elements together with a field-programmable gate array (FPGA) and an external direct current (DC) source 39,40 . Moreover, a series of recent studies demonstrated that metasurfaces composed of Schottky diodes were capable of sensing different waves even at the same frequency in response to the incoming waveform, or more specifically the pulse width [41][42][43][44] . This waveform selectivity provides a new degree of freedom to control electromagnetic waves at a single frequency, which has thus far been exploited to mitigate electromagnetic interference 44 , design antennas [45][46][47] and related microwave devices 48,49 and control communication signals 43,50,51 . Although the electromagnetic response related to the conducting geometry is well known to be improved by existing methods, including the use of equivalent circuit models 52,53 and Jones matrices 54,55 , the relationship between the waveform-selective performance and SPICE parameters of diodes remains unclear, which is important for the design of waveform-selective metasurfaces. For this reason, this study reveals how the waveform selectivity is determined by the SPICE parameters used for waveform-selective metasurfaces. In particular, we present simulation results to show how the SPICE parameters are associated with the operating frequency and power of waveform-selective metasurfaces. Thus, our study provides design guidelines for circuit-based waveformselective metasurfaces to select/fabricate optimal diodes and enhance the waveform-selective performance at the target frequency and power level.

Results
In this study we use the unit cell of the simple waveform-selective metasurface drawn in Fig. 1. In this structure, the unit cell is composed of a square conducting patch, a dielectric substrate (Rogers3003 with loss tangent of 0.0013) and a ground plane (perfect electric conductor: PEC), which resembles unit cells of existing metasurfaces 3 . However, the gap between conducting patches is connected by a set of four Schottky diodes that serve as a full-wave rectifier. Therefore, although the resonant frequency is determined by the same concept as the one that applies to ordinary metasurfaces 3 , the incoming waveform following a sine function is fully rectified to generate another waveform based on the modulus of the sine function. In this case, as seen in the Fourier expansion of the rectified waveform, the incoming frequency component is converted to an infinite set of components, although most of the energy is at zero frequency 41,43,44 . For this reason, the circuit configuration drawn in Fig. 1 enables exploitation of the transient phenomena well known in DC circuits even if the incoming waveform is an alternating current (AC) signal. In the particular case shown in Fig. 1, the energy of a short pulse is temporarily stored in the capacitor and then dissipated with the parallel resistor, resulting in strong absorption. If the pulse width is sufficiently long, however, then the capacitor is fully charged such that electric charges induced by the incident wave cannot enter the internal circuit. Thus, the absorption performance is significantly reduced even at the same frequency. Other types of waveform-selective metasurfaces are further explained in the literature [43][44][45] . In this study, we adopt a co-simulation method available from an ANSYS numerical simulator (Electronics Desktop 2020 R2) to simulate the abovementioned capacitor-based waveform-selective metasurface. Although the co-simulation method is fully explained in the literature 44 , we calculate the absorbing performance of the waveform-selective metasurface for long pulses and short pulses in the frequency and time domains, respectively. For the frequency-domain calculation, we use a small signal analysis (for linear simulations) or harmonic balance analysis (for nonlinear simulations), which readily and effectively estimates the absorptance for long signals. For the time-domain calculation, transient reflected energy is calculated and compared to incident energy, which produces transient absorptance for short pulses. Note that the Schottky diodes used in this study are modelled by referring to the SPICE parameters of a commercial product provided by Avago, specifically the HSMS-286x series (Table 1). However, since we vary these SPICE parameters over a wide range, a similar conclusion is expected to be drawn in designing waveform-selective metasurfaces containing other diodes.
First, we numerically studied how the operating frequency of the waveform-selective metasurface can be scaled to a smaller or larger frequency range. Here, we changed the physical design parameters to build simulation models near 280 MHz, 2.8 GHz, 28 GHz and 280 GHz (see the parameters in Table 2). The substrate thickness was adjusted to be almost consistent with the thickness of existing printed circuit board (PCB, Rogers3003) products. Additionally, the incident pulse width was set to 500 ns, 50 ns, 5 ns and 0.5 ns for the 280-MHz, 2.8-GHz, 28-GHz and 280-GHz models, respectively, to maintain the relative bandwidth of the incident pulse (see the literature reporting that the spectrum spreads if the pulse width is too narrow 56 ). Before evaluating the waveform-selective response, we simplified the waveform-selective metasurface models to replace all the circuit components with a single 120πresistor and investigate the operating frequency by using the abovementioned co-simulation method, more specifically the small signal analysis. In this way, the use of the single resistor was useful for clarifying the relationship between the operating frequency and the physical dimensions of the metasurface models as well as the influence of SPICE parameters used later. The simulated absorptances are plotted in Fig. 2, where the absorptance peak appeared at approximately 280 MHz, 2.8 GHz, 28 GHz and 280 GHz, as expected.
After obtaining these results, we conducted a transient analysis to send sufficiently large signals (0 dBm) to the four models. The absorptance of the waveform-selective metasurface was obtained, as shown in Fig. 3. According to this figure, the 280-MHz model and the 2.8-GHz model showed clear differences between the absorptance for short pulses and that for continuous waves (CWs) due to the abovementioned waveform-selective absorption mechanism. However, the 28-GHz model and the 280-GHz model exhibited almost no difference despite the different waveforms and the presence of the rectifying circuit. This occurred because the diodes contained a Since the reactance of C j approached zero and the diodes no longer rectified induced electric charges, the zero frequency component was severely suppressed in the 280-GHz model. Moreover, the C j used in SPICE parameters was varied, as shown in Fig. 5, where the difference between the pulse absorptance and CW absorptance was found to increase if C j was sufficiently small. For instance, when C j = 1 fF, the 280-GHz model exhibited strong absorptance of approximately 0.9 for a short pulse and limited absorptance of less than 0.2 for a CW. Therefore, reducing C j is very important to distinguish a short pulse from a CW in the high-frequency range. Also, note that this femtofarad-order parasitic capacitance is too small as ordinary commercial products but is realizable as customized integrated circuit chips 57 .
Another point here is that the 280-GHz model had a very small gap of 0.015 mm between conducting patches. This indicates that the loaded circuit components, including diodes, must be smaller than 0.015 mm. From the      58 . To further reduce the chip dimensions, directly fabricating conducting patterns of metasurfaces on chips, namely, not on PCBs, is more realistic 59,60 . Additionally, C j can be reduced to the femtofarad order, as plotted in Fig. 5 59 . However, this does not necessarily ensure that other diode characteristics such as optimal turn-on voltage and current are realized in this power range. Therefore, reaching a balance between C j and other parameters is important to optimize the performance of waveform-selective metasurfaces in the high-frequency range. Next, we investigated how the operating power level is determined by the SPICE parameters. To readily understand the relationship between the operating power level and diode properties, first, we plotted simple IV curves of a single diode biased by a DC source, as drawn in Fig. 6a, where we changed saturation current I S and breakdown voltage V B . According to Fig. 6b, the turn-on voltage was decreased by increasing I S . At the same time, however, increasing I S allowed an increase in the backward current. Additionally, Fig. 6c shows that V B simply determined the voltage to break down the diode. Specifically, the larger V B was, the larger the voltage that could be applied across the diode without current leakage.
Based on the simulation results in Fig. 6, we then varied the I S and V B of the waveform-selective metasurface, as plotted in Fig. 7. Note that in this figure we used the 2.8-GHz model simulated in Fig. 3. Figure 7 shows that by increasing I S , strong absorption was obtained for short pulses with a small power level (e.g., near -10 dBm with I S set to 0.5 nA). At the same time, however, the CW absorptance gradually increased since, as shown in Fig. 6b, the diodes permitted more electric charge to come in from the backward direction as well. In addition, Fig. 7 indicates that by increasing V B , the CW absorptance was suppressed for a wider range of the input power. This figure shows that if the CW absorptance starts increasing, then the input power exceeds the breakdown voltage of the diodes. We also varied other parameters such as forward knee current I K and backward knee current I KR , each of which is related to the amount of current at the turn-on voltage and the breakdown voltage, respectively. Changes in I K and I KR also determine the waveform-selective response, as shown in Fig. 7.
In Fig. 7, we clarified how the operating power range was associated with I S and V B . However, the difference between the pulse absorptance and the CW absorptance was maximized only in a narrow power range (e.g., around 0 dBm in Fig. 7b). This occurred because this dynamic range was constrained by the change in the resistive component of the diodes 61 . More specifically, the entire metasurface impedance including the diode resistance matched the wave impedance of free space to maintain the large absorptance for short pulses in a limited   Fig. 8a. In this case, the resistive component of the diodes still varies with the input power. However, the influence can be almost ignored if R add is moderately large. For this reason, the difference between the pulse absorptance and the CW absorptance was maintained at 0.7 or higher from 0 to 14 dBm when R add was set to 200 , as shown in Fig. 8b. Thus, the dynamic range of the waveform selectivity can be readily broadened by adjusting R add . Moreover, the parasitic series resistive component of the diodes R S relates to the dynamic range of the waveform-selective metasurface to some extent. This is seen in Fig. 9, where R S was changed from 0 to 100 . The absorptance for the short pulse was as large as 0.9 between 0 and 20 dBm when R S was 100 , which was similar to the pulse absorptance with R add = 200 in Fig. 8b. However, the absorptance for the CW was also increased in the same power range because R S was involved in not only the rectification process required for the waveform selectivity but also the intrinsic resonance mechanism of the patch structure (see the structure showing large absorptance in Fig. 2). Therefore, including a resistive component inside the diode bridge as R add is important to increase the absorptance for short pulses but not CWs.
This study numerically unveiled how SPICE parameters of diodes are related to the performance of waveform-selective metasurfaces that are capable of varying electromagnetic response in accordance with the pulse width of incoming waves even at the same frequency. This study only focused on numerical evaluation since measurements add additional complexities to the assessment (e.g. additional parasitic parameters). We note that our simulations were fully based on the co-simulation method of Ansys Electronics Desktop, namely, one of commercial software products widely used and verified by measurements already 41,43,44,47,49,50 . However, the experimental validation is still important to assess the practical performance of waveform-selective metasurfaces with various SPICE parameters, which can be further explored as a future subject. Here we note that the design of the entire circuit configuration is highly customizable. For instance, our study assumed R add as a discrete lumped component, which in practice leads to adding an additional inductive component due to the soldering. However, this parasitic component can be suppressed by integrating R add with diodes as a single packaged chip. In other words, additional parasitic components appearing in measurements highly depend on the circuit design.

Conclusion
In this study, we investigated how the absorption profiles of waveform-selective metasurfaces are related to their SPICE parameters. In particular, we have drawn conclusions about how the SPICE parameters are related to (1) the high-frequency operation, (2) input power requirement and (3) dynamic range of waveform-selective metasurfaces with supporting simulation results. First, we clarified the relationship between the SPICE parameters and the operating frequency. The operating frequencies of conventional metasurfaces that do not include circuit components are well known to be adjusted by simply varying the physical dimensions. In the case of waveform-selective metasurfaces, however, junction capacitance C j was shown to play an important role in the high-frequency range in maintaining the difference between the absorptance for a short pulse and that for a CW. Second, we showed how the operating power level of the waveform-selective metasurfaces is determined by the SPICE parameters. Our simulation results showed that saturation current I S and breakdown voltage V B are important for reducing the operating power level and maintaining small absorptance up to a large input power level, respectively. However, these two parameters were found to not contribute to broadening of the dynamic range of the waveform-selective metasurfaces, which can be improved by introducing an additional resistive component R add into the inside of the diode bridge. Thus, our study is expected to provide design guidelines for circuitbased waveform-selective metasurfaces to select/fabricate optimal diodes and enhance the waveform-selective performance at the target frequency and power level. Our results are usefully exploited to ensure the selectivity based on the pulse duration of the incident wave in a range of potential applications including electromagnetic interference, wireless power transfer, antenna design, wireless communications, and sensing.